Programmed buckling by controlled lateral swelling in a thin elastic sheet 
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Recent experiments have imposed controlled swelling patterns on thin polymer films, which subse- 
quently buckle into three-dimensional shapes. We develop a solution to the design problem suggested 
by such systems, namely if and how one can generate particular three-dimensional shapes from thin 
elastic sheets by mere imposition of a two-dimensional pattern of locally isotropic growth. Not every 
shape is possible. Several types of obstruction can arise, some of which depend on the sheet thick- 
ness. We provide some examples using the axisymmetric form of the problem, which is analytically 
tractable. 



I. INTRODUCTION 

The inhomogeneous growth of thin elastic sheets is 
emerging as a powerful method for the design of three- 
dimensional structures from two-dimensional templates 
0, ■ Most of the associated theory has focused on pre- 
dicting the buckled shape of a sheet that results from a 
given imposed pattern of growth or swelling 3-9] . This 
work has led to a number of insights, and continues 
to present challenges for theorists and experimentalists 
alike. For most applications, however, one knows the de- 
sired surface but not the growth pattern that generates 
it. In this paper, we pose this reverse buckling problem 
and solve its axisymmetric form. 

The reverse problem has received little attention thus 
far, and one might mistakenly assume that the solution is 
trivial. Certainly, it is an elementary exercise in differen- 
tial geometry to determine the unique metric associated 
with a surface. With recourse to through-thickness varia- 
tions in material properties, extrinsic curvatures may also 
be programmed to select a unique shape. However, real 
material systems may be more limited. Thus, we wish to 
explore the question of which shapes can be made with 
thin elastic sheets when only the two-dimensional mid- 
surface metric can be prescribed. In a practical sense, 
these are the shapes programmable by encoding a single, 
spatially-dependent scalar property into a sheet, namely 
an isotropic swelling ratio [lOj. 

How could programming the metric of a shape fail to 
reproduce the shape? A metric may have families of 
shapes from which to choose its buckled configuration. 
One sees this possibility immediately using a flat metric; 
a continuum of cylindrical, truncated conical, and other 
developable immersions exist, but a piece of the plane has 
lower bending energy than all. Beyond this, the actual 
metric realized by the sheet often differs from the pre- 
scribed metric because accommodating a finite-thickness 
bending energy may induce an in-plane strain through 
geometric compatibility conditions. Finally, there is no 
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guarantee that a given metric can satisfy all of the bound- 
ary conditions at the sheet edges. Such a sheet presum- 
ably forms a boundary layer which, while vanishing in 
the zero thickness limit [8|, may have nontrivial effects 
at finite thickness. 

Thus, we wish to know how to prescribe a metric on a 
sheet of given thickness to produce a desired shape, and 
what limits there are to the shapes that can be prescribed 
exactly, up to and including their boundaries. We will 
make some progress towards answering these questions 
in what follows. 



II. EQUILIBRIUM SHAPE EQUATIONS 

Consider a flat sheet and multiply all of its metric 
components by a nonuniform scalar field that repre- 
sents the degree of local, isotropic swelling. The sheet's 
preferred extrinsic curvature remains zero locally every- 
where, but we have prescribed a new "target" [11| metric 
a a p (Greek indices range from 1 to 2) on its midsurface. 
The actual metric a a p adopted by the sheet's midsurface 
will, in general, not agree with a a $ for sheets of finite 
thickness. However, when designing a shape, it is this 
realized metric that is known. We suggest that the nat- 
ural formulation for our design problem is to begin with 
the realized metric a a p of a desired shape in its most 
convenient coordinate system, and then compute the ap- 
propriate target metric a Qj g. Finally, we must return to 
a coordinate system convenient with respect to the flat 
configuration. What emerges is not only a swelling fac- 
tor f2, if one exists, but a set of boundaries that may 
be cut into the unswelled sheet to produce the desired 
boundaries of the final shape. 

In the appendices, we derive our equilibrium equa- 
tions from a three-dimensional elastic energy, following 
the procedure of Efrati and co-workers pj. We have 
adapted their approach to our purpose by using the re- 
alized metric, a a p, rather than the target metric a a p, 
to raise and lower indices as well as to construct a co- 
variant derivative, a measure of integration, and an elas- 
tic tensor. As long as the three-dimensional strains are 
sufficiently small, the difference in our three-dimensional 
energies is higher than quadratic order in these strains. 
The result of this change of viewpoint will be a set of 
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equations involving the two curvature invariants of the 
surface, with an immediate resemblance to well-known 
equations describing fluid membranes fl^ . 

Consider a thin sheet of thickness t with midsurface im- 
mersion X and unit normal N. The metric and curvature 
tensors of the surface X are the first and second funda- 
mental forms a a p = <9 Q X ■ 9^X and b a p = d a dpX. ■ N = 
— 9 Q N • <9^X. Two independent invariants of the curva- 
ture tensor may be represented using geometric quanti- 
ties, namely twice the mean and Gaussian curvatures: 
2H = 6° and 2K = b%b% - fo^fog. The prescribed met- 
ric a a p corresponds to the metric of an "unstretched" 
midsurface. The deviation of the realized metric from 
the prescribed metric defines twice the in-plane strain: 
2e Q( 3 = a a p — a a j3. 

In the limit of small thickness, we may write a two- 
dimensional energy (Appendix |A|) : 
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is a two-dimensional elastic tensor incorporating an 
isotropic Young's modulus and Poisson's ratio Y and v. 
The first and last terms in this energy are the stretch- 
ing and bending terms of Efrati et al. @|. We view t\\b\\ 
as our small parameter, but expect ||e|| to be small near 
equilibrium. Accordingly, we retain three coupled strain- 
curvature terms of order i 3 ||o|| 2 ||e| dropped in [7], but 
not terms of order i 3 ||o|] 2 ||e|| 2 . One may ask why such 
coupled terms should be retained, if when s is small they 
contribute to the energy at higher order. We will see that 
the corresponding terms in the Euler-Lagrange equations 
of equilibrium are independent of ||e||, and of the same 
order as terms arising from variation of the i 3 ||6|| 2 bend- 
ing term. Similarly, the stretching term of order i||£|| 2 
produces terms of order t||e|| in the equations. 

Variation of the immersion of the midsurface (Ap- 
pendix [B]) leads to equations of equilibrium 

2B [V a V a H + 2H (H 2 - K)} - s a0 b afi = 0, (3) 

V o a a/J =0, (4) 

free smooth boundary conditions [13] 

2n a V a H + (1 - v^Vib^nJ 

n p [B (2H 2 - (1 - u)K) a al3 + s 

n a nfi [2vHa afi + (1 - v)b a ? 

and a corner jump condition for free piecewise-smooth 
boundaries 
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Here, V Q is a covariant derivative constructed from the 
realized surface metric, n and I are surface tangents nor- 
mal and tangent, respectively, to the boundary, [J de- 
notes a jump in the enclosed quantities, and B = 12 (Y— 1^2) 
is a bending modulus. We have neglected terms of orders 
t 3 ||6|| 2 ||e|| and i||e|| 2 , with the implicit assumption that 
derivatives do not affect order. Finally, the "effective" 
stress tensor is given by 

ga p = t A a ^ b e 1& + (HA a ^ s b l5 - b a K A K ^ s b l8 ) 



-A a ^ s b;b K 



(9) 



It should be apparent from this expression that an un- 
stretched midsurface, that is, one free of in-plane strain, 
does not imply an unstressed finite-thickness sheet. This 
extrinsic contribution to the stress is the only result of 
our retention of coupled strain-curvature terms in the 
energy. The tensors tA aM& £ lS and ^A a ^ 5 b lS are the 
stress and moment tensors of Efrati et al. Q; after ex- 
plicitly raising indices, the latter becomes the bracketed 
term in the torque boundary condition (J?) and corner 
condition ©. The final term in the definition of the ef- 
fective stress is an analogous application of the elastic 
tensor to the "third fundamental form" c^s = b*b K s- 

This stress plays the same role as that of the Lagrange 
multipliers of Guven & Miiller 14|, whose equilibrium 
equations for paper coincide with ours when K = 0. Ac- 
tually, the definition of "stress" is rather malleable. Our 
definition's inclusion of the extrinsic terms that manifest 
in the divergence-free quantity in Q seems natural, es- 
pecially as they arise from variation of the energy with 
respect to the in-plane strain. Operationally, any scalar 
T that does not vary, or varies such that ST — T Q ^<5a a/ 3, 
will produce only terms that may be tucked into s Q/3 . We 
note also the absorption of gravitational forces into the 
Lagrange multipliers in fl4| , reminiscent of the definition 
of "dynamic pressure" in problems involving isochoric flu- 
ids. For more ambiguities, see [l5| and the discussion of 
"null stresses" in jl6j |. 

If these equations are to be solved for the six terms a a p 
and b a fj, rather than an explicit immersion X, they must 
be supplemented by the Peterson-Mainardi-Codazzi and 
Gauss equations: 
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.(ii) 



The rjg 7 are the usual Christoffcl symbols. These auxil- 
iary equations are automatically satisfied by any immer- 
sion X or X + <5X. 

Many shapes cannot satisfy the boundary conditions 
([5][7|). For example, the normal force and torque bound- 
ary conditions ([5]) and ([7]) are incompatible for minimal 
surfaces of the helicoid-catenoid family. Such shapes re- 
quire either a boundary layer or applied boundary forces 
and torques. 
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A consequence of the corner condition (|8j may be ob- 
served by bending two adjacent sides of a piece of paper 
towards each other. Curvature must vanish "across" the 
sheet at the corner, so the tip remains flat. 



III. MAKING SHAPES 

When the prescribed metric is given and one is solv- 
ing "forwards" for the shape, one must solve the equilib- 
rium equations (J3J4J) and geometric integrability condi- 
tions (jlOtilll) for the six components a a /3 and b a p . In the 
reverse problem, we choose a shape X that satisfies two 
boundary conditions, ([5J) and ([7]); integrability is auto- 
matically satisfied. After solving the equilibrium equa- 
tions ([3]|3]), along with boundary condition ©, for the 
components of the stress tensor s Q/3 , we recover via (|9]) 
the target metric a a p in whatever buckled coordinate sys- 
tem we chose for our initial convenience on X. Finally, we 
must determine a coordinate transformation back into an 
appropriate laboratory frame for assigning the swelling 
factor to the unbuckled sheet [17| . 

Though many shapes do not satisfy all of the boundary 
conditions, in principle only a boundary layer is needed 
to balance the normal force (|5|) and torque ([7]) condi- 
tions. Note that this layer may be incorporated into the 
prescribed swelling factor, so need not share the char- 
acteristic width of spontaneously formed layers @. The 
tangential force condition (J6]) is more involved. In gen- 
eral, the integration constants of the first-order equation 
(U]) may be insufficient to balance these in-plane forces, 
which may require global changes in the metric. Below, 
we will explore how these boundary conditions affect the 
construction of axisymmetric shapes, for which we solve 
the equations of equilibrium analytically. 



and one algebraic equation 
d 2 p 
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We can integrate equation (|13l) to yield 



(15) 
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dyg(y)p{y) yi - [d y p(y)Y 



. (16) 

where C is an integration constant, and Ub > u lies on 
one boundary of the sheet. Given the stress tensor, Eq. 

is now an algebraic equation for the strain tensor and, 
thus, the prescribed metric a a p. 

Our coordinates (u, v) are natural for the buckled ob- 
ject, but not for the laboratory. We must perform a 
change of variables to a coordinate system convenient for 
programming an isotropic swelling factor SI (r). A nat- 
ural choice for axisymmetric shapes is to use cylindrical 
polar coordinates (r, 0) and identify v with 9, so that 
the metric becomes Vt{r)(dr 2 + r 2 d9 2 ). The coordinate 
transformation u{r) is determined by the solution to the 
differential equation 



[d r u(r)Y 



and the swelling factor by 



a V v[u(r)} 
r 2 a uu [u(r)] ' 



fi(r) 



a vv [u(r)j 



(17) 



(18) 



A. The axisymmetric solution 



B. The axisymmetric boundary conditions 



Consider an axisymmetric shape parametrized by arc 
length along a meridional geodesic, so that: 

X(u,v) = ^p{u) cos v, p(u) sin v, J dy \J\ - [d y p{y)} 

(12) 

Thus, a a pdx a dx^ — du 2 + p 2 dv 2 and b a pdx a dx^ — 

- d 2 u pi^/\-{d u p) 2 \du 2 + p^\-{d uP ) 2 dv 2 , with 

mean curvature H = [1 — d u (pd u p)] / 2py/l — {d u p) 2 
and Gaussian curvature K — —d 2 p/ p. 

We are free to specify the function p, as long as 
(d u p) 2 < 1. If we assume a diagonalized, axisymmetric 
target metric, the equilibrium equations reduce to one 
differential equation 



+ s uu d u In py/i - (d uP ) 2 + g(u) = , (13) 



Our construction relies on our ability to find an X 
that satisfies the free boundary conditions. Two of these, 
([5]) and (|7|). are simply two conditions on p{u) on the 
boundaries. On the boundary Ub, the first of these is 



p"{ub) = v 



1 - P l2 {u b ) 
P(ub) 



(19) 



which implies K(uY) = — p h '(ub) / p(ub) < 0. The second 



p"'{u b ) = + ^ p'{ub)K(ub) . (20) 

These conditions can be easily satisfied using an arbi- 
trarily narrow region near the boundary u b - If there is 
another boundary u a < Ub, the same considerations ap- 
ply there. 
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(a) Swelling factor 




(b) Ziggurat 

FIG. 1. (Color online) (a) Swelling factor Q(r) that swells 
a disk into the shape in (b) for thicknesses t = 1/100 (solid 
black), t = 1/200 (dashed black), and t = 1/500 (solid grey). 
The outer radii needed are w 0.9, « 1.2 and « 1.3, respec- 
tively. Lengths are in units of the radial arc length of the final 
shape. 

We can also satisfy the boundary condition (|6|) at u b 
by choosing the integration constant 

C = -B [2H(u b f - (1 - v)K(u b )\ p{u b )^l-p'{u b f. 

(21) 

If the surface has only one boundary, this procedure is 
sufficient. However, with two boundaries we must also 
satisfy an integral constraint, 

rut 

/ du g{u)p{u) Vl - P'W 2 = (22) 

-C-B [2H(u a ) 2 - (1 - v)K{u a )] p{u a )^\- p'{u a f . 

This is a nonlocal constraint as it involves both bound- 
aries, located at u a and u b . This global balance may 
require a global change in p(u) to accommodate. Our 
procedure, by no means unique, for this accommodation 
is shown below in the example of the asymmetric annular 
sheet. 



C. Examples 

We start our survey of examples with a topological 
disk. The swelled shape is shown in Fig. [T^b) and the 



swelling factor used to produce it is shown in Fig. Ufa) 
for three different sheet thicknesses. This figure was pro- 
duced using a metric of the form 

p(u) = u + A lU 4 e~ 10 [e 10u - l] (23) 

0.3m 5 

~ 0.01 + u 4 /[sin(77ra)/(77r) + 0.5] ' 

where A\ and A 2 were chosen to satisfy the normal force 
and torque boundary conditions at u = u b = 1. We find 
Ax sa -0.829 and A 2 « -0.723. Despite its seeming 
absurdity Eq. (|23p underscores the flexibility we have 
in choosing a metric. Moreover, there is some method 
to our choice. Since stresses cannot diverge, we require 
that p(u) asymptotically flatten at the center. It is a 
straightforward calculation, by expanding the stresses in 
a power series in u, to show that this requires p(0) = 
0, p'(0) = 1, p"(0) = 0, p"'(0) < 0, and p""(0) = 0. 
Following the procedure described in the previous two 
sections, the swelling factor is easily obtained. 

The process of choosing coefficients for Eq. (|2"3")l reveals 
some of the potential pitfalls of designing a shape. Some 
Ai and A 2 that satisfy the boundary conditions require 
|p'(u)| > 1 at one or more places within the sheet. More- 
over, one could find that the resulting prescribed metric 
is not positive definite everywhere within sheets that are 
too thick. 

We now consider a pair of topological annuli, shapes 
with two free boundaries. Satisfying conditions on both 
boundaries is simple for a surface symmetric about a fixed 
u. Choosing C to satisfy one boundary, we automatically 
satisfy the other. An example is shown in Fig. [5] The 
metric we use is 

p(u) = A ie - 10 [e 10u - 1] + A 2 (l - u)e- w [e Wu - l] 

+Bl [e~ Wu - e- 10 ] + B 2 u [e- Wu - e- 10 ] (24) 

i • IK \ 

+1 H sin(o7ra). 

57T 

We find Ai = Bi w 0.205 and A 2 = B 2 w 1.013 in 
order to satisfy the normal force and torque boundary 
conditions at u = u a = and u = u b = 1. 

An asymmetric annulus is significantly more compli- 
cated. It is no longer sufficient to simply choose a met- 
ric appropriately on the boundaries, because the integral 
constraint (f2"2"| depends on the value of p(u) throughout 
the sheet. After satisfying one boundary, it will generally 
be impossible to satisfy the other without modifying the 
metric. For the example in Fig. [31 we use 

p(u) = A ie ~ w [e Wu - 1] + A 2 (l - u)e~ w [e Wu - l] 

+Bl [ e - Wu - e- 10 ] + B 2 u [e- Wu - e" 10 ] (25) 

_^ + .5+-e- 32 ("- 3 / 4 ) 2 , 
16 

adjusting A\, A 2 , B\ and B 2 according to the normal 
force and torque boundary conditions at u = u a = 
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(a) Swelling factor 




(b) Sheave 

FIG. 2. (Color online) (a) Swelling factor f2(r) that swells an 
annulus with inner radius r — 0.5 and outer radius r ~ 1.28 
into the shape in (b) for thicknesses t — 1/100 (solid black) 
and t = 1/500 (dashed black). Lengths are in units of the 
radial arc length of the final shape. 

and u = Ub = 1, choosing C to satisfy the tangential 
force boundary condition at u = it& = 1 and the param- 
eter rj to satisfy the tangential force boundary condition 
at u = u a = 0. We find A x m -0.044, A 2 « 0.029, 
B\ w —0.18 and B 2 ~ 0.10. Although there is some 
weak thickness dependence in r], we find ?/ w —0.051 for 
thicknesses from t — 1/20 to t — 1/500. There is also a 
weaker dependence on thickness for fl(r) in this example, 
presumably because this swelled shape has less curvature 
than the others. 

Our final example is the disk shown in Fig. m with 

p(u) a (6.22 x 10" 11 ) e" 10 (e Wu - l) u 4 

+ (-2.1x 10- n )e- 10 (e 10u -l) u 5 

+ 1 [(5.966 - 7.02w)erf(7.07u - 6.01) 

8 

+ (u[1.2u - 2.16] + 0.9855) erf(6.67w - 6) (26) 

0.13u 3 (3 3Q x ^ u4 L2ua 

125w 3 + l v ; 
-1.18« + 6.95 + e^ 44 - 44 ^°- 9 ) 2 (0.1 U - 0.09) 

-o.see- 50 ^- - 85 ) 2 " 

This "drum" requires a disk of radius 60 to produce, 
though it has a total center-to-edge arc length of only 
2.75. This requires significant shrinking, up to a local 
swelling factor of ss 5 x 10~ 5 . This, of course, is another 
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(a) Swelling factor 




(b) Compression fitting 



FIG. 3. (Color online) (a) Swelling factor fi(r) that swells an 
annulus with inner radius r — 0.3 and outer radius r ~ 1.81 
into the shape in (b) for thicknesses t — 1/20 (dashed black) 
and t — 1/100 (solid black). Lengths are in units of the radial 
arc length of the final shape. 

obstruction to swelling a shape: the required swelling 
factor may be beyond the capabilities of any existing ex- 
perimental system. So one can swell the shape of a drum 
in principle, though perhaps not currently in practice. 
We note, however, that conformal transformations of the 
prescribed metric into some atypical coordinate system 
may be a way to improve the range of required swelling 
factors. 



IV. DISCUSSION AND CONCLUSIONS 

We have presented equations that govern the design of 
shapes by isotropic growth of a thin elastic sheet, and 
solved them analytically for axisymmetric cases. There 
are two relevant cases: disk-like sheets with one bound- 
ary, and annular sheets with two boundaries. For disk- 
like sheets, the boundary conditions can be satisfied lo- 
cally by choosing appropriate gradients in the metric near 
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(a) Swelling factor 
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FIG. 4. (Color online) (a) Swelling factor Q(r) that swells 
a disk with outer radius r ~ 60 into a "drum" with a total 
radial arc length of 2.75 units, for thicknesses t — 1/100 (solid 
black), t = 1/200 (dashed black), and t = 1/500 (solid grey). 
We show f2(r) only up to r = 1.5, beyond which it simply 
approaches zero. 



the edges. This implies, among other things, that a neg- 
ative Gaussian curvature lip must appear. For a generic 
annular sheet, not only must we satisfy the local bound- 
ary conditions at two boundaries, but a difficult nonlocal 
condition resulting from in-plane force balance. We have 
found that an additional term in p{u) linear in u with 
adjustable coefficient can be used to satisfy this bound- 
ary condition without changing the surface dramatically. 
Once the assumption of axisymmetry is lifted, it is not 
at all clear what would be required to satisfy this last 
boundary condition. 

Several obstacles to swelling a shape may arise. A cho- 
sen functional form for p(u) may, after boundary con- 
ditions are applied, fail to satisfy |p'(it)| < 1 at one 
or more places on the sheet. This implies that an ax- 
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FIG. 5. Nondimensionalized azimuthal stress, s vv /(tY), in 
the swelled surface of Fig. [3] with t = 1/20. 



isymmetric shape of the desired form cannot satisfy the 
boundary conditions. For a sufficiently thick sheet, the 
required prescribed metric may fail to remain positive 
definite at one or more places. Thus, the shape would 
not be swellable, at least not by an axisymmetric tar- 
get metric, at the desired thickness. This is not a sur- 
prising result, as sufficiently thick sheets will not buckle 
at all in response to an 0(1) inhomogeneity in swelling. 
Yet another problem- one that we did not encounter- 
may arise because our method of search for an isotropic 
swelling factor rests on an implicit assumption about the 
existence of a global conformal coordinate system. Such 
a coordinate system is only guaranteed to exist locally 
[l8| . Failure to find a coordinate system of this type glob- 
ally suggests, again, that a particular shape may not be 
swellable. Finally, we note that we have not investigated 
the stability of the generated surfaces, only whether ax- 
isymmetric extrema exist. The azimuthal stress s vv (Fig. 
[5} in the surface from Fig.[3]oscillates between tensile and 
compressive, which suggests the possibility of a wrinkling 
instability in the compressed regions of sufficiently thin 
sheets. 
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Appendix A: Derivation of the 2D Energy 

Our procedure mirrors that of Efrati et at. Q, but 
involves some additions. In what follows, Latin indices 
run from 1 to 3, Greek indices from 1 to 2. 

We begin with a three-dimensional body with elastic 



7 



energy 



E 



(Al) 



= \j y/9 [A(Tr(e)) 2 +2 A ,Tr( e 2 )] , 

using the following definitions for the elastic tensor and 
strain: 



2e, : 



Y 



l + v 



9ij 9ij 



Here Y and v are an isotropic Young's modulus and Pois- 
son's ratio, and g~ij is a target metric. 

If our body is a thin plate or shell, we may view it as 
a stack of surfaces, express quantities on each surface as 
expansions around the middle surface [l9j . and integrate 
along the thin dimension z to obtain a two-dimensional 
energy. A plate is merely a special case of a shell, in 
which the target metric is independent of z. 

The position vector in the body is expressed using the 
midsurface immersion and unit normal 



R = X + zN , 

so the metric is block-diagonal 



= a-R-<9,R = 



(9af3) 2 x2 O 

O i 



(A2) 



(A3) 



and tangents to surfaces in the stack are 
d a R = d a X+zd a N = (££ - zb'l) dpX = i&dpX , (A4) 

where we have defined a tensor 7r^. This tensor has an 
inverse pp" such that 



d a R = pJJ0"X 



and hence 



d a R ■ dfjK = g a0 = n^pa^ . 
d a R ■ d R = 5^ - - Q - 7 
d a R ■ d p R = g a0 - 
With these definitions, we find that: 



jpn. = Op = p^'p , 

P?pf a 74 ■ 



+ zbi + zWbl + 0(z 3 ) 



and 



^={l-2zH + z 2 K) s fa. 



(A5) 

(A6) 
(A7) 
(A8) 

(A9) 
(A10) 



Movingfrom geometric to physical quantities, we apply 
the first [3] Kirchhoff-Love assumption A lZkl eki — to 
express our elastic energy as 



A 11 eijeki 



Y 



l + i/ 



e2(z)e P p(z) + e a p(z)ei(z) 



A a ^ s (z)e a p(z)e lS (z). 



(All) 



The second [tJ Kirchhoff-Love assumption = gives 
us a block-diagonal form for the target metric: 



'(a Q/3 ) 2x2 O 

9a/3 - J q 1 



(A12) 



Finally, the expansion of the elastic tensor takes the form 



k,l,m,n—0 

= A afilS + 2z (A K ^ 5 b a K + A aPKS bl) + 0{z 2 ) , 

(A13) 

where we have defined (b k )" as the contracted product 
of k 6-components bpb@ . . . b]. and used the definition ((2|) 
from the main text. 

Rolling all of this together, and recalling the definition 
2e a /3 = a a p — a a p, our elastic energy (IA1|1 is 



1 r - r 1 

E=- y/a dz (1 - 2zH + z 2 K) (A14) 

2 Js J-t/2 

x [A aMS + 2z (A KM5 b a K + A af}KS bl) + 0{z 2 )} 
x [sape-fS - 2ze a pb lS + z 2 (b a pb jS + e a pb^b K& ) 
+0(z 3 )] , 
which evaluates to (JlJ. 

Appendix B: Variation of the 2D Energy 

We consider a variation of the midsurface configuration 
SX = Su a d a X + 5(N , (Bl) 
under which the forms vary as follows [20|: 

5a aj3 = W a Sup + VpSua - 2b a p6(, (B2) 
5b a p = {\7p5v? - b}SC)b aj + V Q {b 01 5ut) + V Q V^C . 

Straightforward application of these expressions provides 
the useful variations 

SA a(>w = _ (^7/* «« + A aMx a^ K ) 5a KX , (B3) 

5y^= (V a Su a -2H5()y/^, (B4) 



5H = -V Q V Q <SC 



(B5) 



(2H 2 -K) 5( + 5u a V a H. 



while a laborious calculation making use of the Peterson- 
Mainardi-Codazzi relations (fTU)) yields: 

6 (y/aK) = V Q (Vo [K6u a + {2Ha af} - b afi ) V P S(] ) . 

(B6) 

That this last expression is a divergence, and hence only 
relevant on the boundaries, should be expected from the 
Gauss-Bonnet theorem [2 ll . 



8 



The equations (|3]__|) and boundary conditions ([5][8]) fol- 
low from stationarity of the energy (TTJ) with respect to 
arbitrary 5u a and S( in the bulk and on the boundaries, 
boundary-normal tangent derivatives d n 5C on the bound- 



aries, and boundary-tangent tangent derivatives diSC at 
corners, neglecting terms of orders i 3 ||6|| 2 ||e|| and i||e|| 2 
after the variation. 
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